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The behaviour of two-dimensional patchy particles with 5 and 7 regularly- arranged patches is investigated by 
computer simulation. For higher pressures and wider patch widths, hexagonal crystals have the lowest en- 
thalpy, whereas at lower pressures and for narrower patches, lower-density crystals with five nearest neighbours 
and that are based on the (3^,4,3,4) tiling of squares and triangles become lower in enthalpy. Interestingly, in 
regions of parameter space near to that where the hexagonal crystals become stable, quasicrystalline structures 
with dodecagonal symmetry form on cooling from high temperature. These quasicrystals can be considered 
as tilings of squares and triangles, and are probably stabilized by the large configurational entropy associated 
with all the different possible such tilings. The potential for experimentally realizing such structures using 
DNA multi-arm motifs are discussed. 
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I. INTRODUCTION 

Since the formation of quasicrystals were first reported 
in 1984 by Shechtman et al. for an Al-Mn alloy, ^ many 
systems have been to found to exhibit a quasicrystalline 
phase. Most of these are binary or ternary metallic alloys 
(but never a pure metal). The quest to find quasicrystals 
beyond alloys has led to an increasing number of exam- 
ples, mostly in the field of soft condensed mat ter.!^ Exam- 
ples include dendrimers,'^star polymers,^ micelles^ and 
binary mixtures of nanoparticles.^ All of these examples, 
except for one of the micellar systems, El have dodecago- 
nal symmetry and are often found in regions of parameter 
space close to where crystalline approximants, such as the 
Frank-Kasper <j phase, are observed.*^^^ 

A remaining target is to find a colloidal system that 
can self-assemble into a quasicrystalline structure in the 
absence of an external field (two-dimensional colloidal 
quasicrystals can be induced to form using quasiperi- 
odic light fieldJPKi^. One approach might be to use a 
binary or ternary colloidal mixture, but although com- 
plex crystal structures have been reported for binary 
mixtures ,'13'^ as yet no quasicrystals have been ob- 
served. Another approach might be to use colloids with 
anisotropic "patchy" interactions,^^ where the positions 
of the patches could be used to control the preferred lo- 
cal geometry and hence infiuence the global structure 
formed. Indeed, much progress has been made in de- 
veloping methods to synthesize such types of colloidal 
particles. Furthermore, the first experiments on the 
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novel structures that such patchy interactions can en- 
able the systems to adopt are beginning to appear,!^ as 

well as bei ng sys tematically explored through computer 
simulations.l^^KlII 

Quasicrystals have also been found to form for a vari- 
ety of model potentials in computer simulations. Inter- 
estingly, these are not restricted to mixtures, but can 
also occur for one-component systems. ^^ ^^ Examples in- 
clude isotropic pair potentials with both a maximum and 
a minimum in the potential'^^"^^" micellar models, hard 
tetrahedra^^ and triangular bipyramids,^^ and wateif^ 
and silicon^ bilayers. 

Here, we wish to examine in detail the behaviour of 
two-dimensional patchy particles with 5 and 7 patches 
as possible quasicrystal-forming systems. The reason 
we choose this system is that in a recent study of two- 
dimensional disks with regularly arranged patche s, in- 
triguing behaviour was seen for the 5-patch system .^^^^ 
Hard disks naturally form a hexagonal crystal, and six 
patches reinforce this tendency. For four patches, there 
is a competition between a square crystal, which is en- 
ergetically stabilized by the patchy interactions, and a 
hexagonal crystal, which is stabilized by its higher den- 
sity. However, the situation for particles with 5 regularly 
arranged patches is more complex, since the 5-fold sym- 
metry of the particles is incompatible with crystalline 
order — there is no crystal in which all the patches can 
point directly at those on neighbouring particles. 

In Ref . l35l we found that on cooling (at low pressure), 
crystallization was not observed, because of the frustra- 
tion introduced by the geometry of the patches. The 
configurations generated did though show certain com- 
mon local motifs (those shown in Fig. [T]) and could be 
considered as tilings of squares and triangles. Although 
there was no overall crystalline order, there was some evi- 
dence suggestive of longer-range orientational order; how- 
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FIG. 1. The three local environments that are mainly present 
in the simulated structures. 



potential is recovered. For computational efficiency the 
potential is truncated and shifted at r = Sctlj, and the 
crossover distance in Eq. [T]is adjusted so that it still oc- 
curs when the potential is zero. 

As wel l as studying crystalization in two 
dimensions P^ l ^^ l ^ -^ this patchy particle model has 
also been previously used to study crystallization 
in three dimensions pSHSlI ^nd the self-assembly of 
monodisperse shells .l^^^^ 



ever, because of the relatively small systems sizes consid- 
ered this was not pursued further at the time. Given 
the known tendency of random square-triangle tilings to 
form dodecagonal quasicrystals,^^ we now investigate in 
detail the suggestion made in Ref. i35j that this system 
may form quasicrystals. 

II. METHODS 
A. Potential 



B. Simulation 

Systems of particles were simulated using standard 
Metropolis Monte Carlo (MC) in the NPT ensemble. Pe- 
riodic boundary conditions were applied, and the simu- 
lation box was constrained to be square. The MC moves 
were single-particle rotational and translational moves, 
as well as "volume" moves in which the area of the box 
was scaled, the latter allowing a constant pressure en- 
semble to be simulated. The number of particles N was 
always 2500. 



The model consists of spherical particles patterned 
with attractive patches. They are described by a modi- 
fied Lennard- Jones potential, in which the repulsive part 
of the potential is isotropic, but the attractive part is 
anisotropic and depends on the alignment of patches on 
interacting particles. Specifically, the potential is de- 
scribed by 



"•^^ \Vi^j{rij)Vs,ng{rij,ni,nj) r > (JLJ, 

(1) 

where Vlj, the Lennard- Jones potential, is given by 



VLj(r) =4e 



(2) 



The minimum of this potential is at Teq — 2 / o"lj- 
is an angular modulation factor that depends on the ori- 
entations of the patches on the two interacting particles 
with respect to the interparticle vector. Specifically, 



where 



Gij{rij,ni) exp 



(4) 



cTpw is a measure of the angular width of the patches, 
6kij is the angle between patch vector k on particle i and 
the interparticle vector r^j, and /Cmin is the patch that 
minimizes the magnitude of this angle. Hence, only the 
patches on each particle that are closest to the interpar- 
ticle vector interact with each other, and T4ng = 1 if the 
patches point directly at each other. One feature of this 
potential is that as (jpw oo the isotropic Lennard- Jones 



C. Structural analysis 

We employ a number of approaches to analyse the 
structures that the systems of patchy particles adopt. In 
the condensed state, we find that the 5- and 7-patch par- 
ticles virtually always adopt one of three local environ- 
ments. These are illustrated in Fig. [l] There are two pos- 
sible five-coordinate environments, which are based upon 
two different ways of locally packing squares and trian- 
gles. We refer to them as the a and H environment s by 
analogy to the Frank-Kasper phases of these names. 
These a and H Frank Kasper phases are crystalline ap- 
proximants to dodecagonal quasicrystals and can be con- 
sidered as square-triangle tilings in two of their three 
dimensions. The tilings containing only these local en- 
vironments are more formally denoted by (3^,4,3,4) (a) 
and (3^,4^) (JY), where this nomenclature refers to the 
sequence of polygons around a vertex. The third lo- 
cal environment is the six-coordinate hexagonal environ- 
ment, which represents the densest local packing, and 
will be denoted Z, again by analogy to a Frank-Kasper 
phase. 

To identify these environments we employ a common 
neighbour analysis. Neighbours are defined as all parti- 
cles within a certain cutoff radius Tc from a given par- 
ticle. Then for each pair of neighbours the number of 
neighbours common to both is determined. Each local 
environment is characterized by a unique signature in 
terms of the number of neighbours with which the cen- 
tral particle shares a certain number of neighbours. For 
example, a particle in a <j local environment has a total of 
five neighbours, with one of which it has two neighbours 
in common and with four of which it has one neighbour 
in common. Hence, this environment is denoted by the 
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common neighbour signature {21111}. In a similar way, 
the signatures for the H and Z local environments are 
{22110} and {222222}, respectively. Particles that are 
not in any of these three local environments are labelled 
'undefined' (/7). The cutoff radius {vc ^ 1.38 (Jlj) was 
chosen such that the fraction of [/-particles was mini- 
mized. 

To probe the global order of the configurations gen- 
erated, the associated diffraction patterns were calcu- 
lated. Quasicrystalline configurations are characterized 
by diffraction patterns with non-crystallographic rota- 
tional symmetries. We calculated the diffraction patterns 
by evaluating the real part of the interference function: 



N-l N-1 



(5) 
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where q is the wavevector and the position of particle 
i. Unless stated otherwise, the resolution we used for the 
plots of the diffraction patterns is /S^q^ = Ag^ = 0.15 cr^j . 

When mapping out the behaviour of our system as 
a function of the parameter space it will be useful to to 
measure the degree of twelvefold symmetry of the diffrac- 
tion pattern. To achieve this, we take a Fourier transform 
around a ring of the interference pattern at a g value cor- 
responding to the first diffraction peaks. Specifically, we 
evaluate 



^(^) = ^ S{qist, (t)i) exp {-27i\v(l)i 



(6) 



where the sum is over the n sampled points of S{qisti 
A twelvefold symmetric diffraction pattern will have high 
values at v = 12 and multiples thereof. Similarly, a six- 
fold symmetric diffraction pattern, as would be expected 
for a hexagonal crystal, will have high values at = 6 and 
multiples thereof (including v = 12). Therefore, we used 
\F{iy = 12)1 — \F{v = 6)1 as a measure of the twelve- 
fold character of the diffraction pattern. We note that 
the peaks at multiples of the lowest frequency will have 
a lower amplitude than the main lowest frequency peak, 
where this decay in amplitude is stronger if the diffrac- 
tion peaks are more diffuse. 



III. RESULTS 

We will first consider the 5-patch particles in detail. 
The behaviour of the 7-patch particles is fairly similar 
and so we will consider this system more briefiy in Section 

niBi 



A. 5-patch particles 



1. Low-enthalpy structures 



The most stable phase at zero temperature is simply 
that with the lowest enthalpy. Therefore, we minimized 
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FIG. 2. Zero-temperature phase diagram showing the depen- 
dence of the lowest- enthalpy structure on pressure and patch 
width. 



the enthalpy for a variety of candidate crystal structures 
under different conditions and potential parameters in 
order to obtain Figure [2j which shows how the lowest en- 
thalpy crystal structure depends on pressure and patch 
width. This figure is slightly different from that which 
appeared in Ref. [35l because a slightly lower energy crys- 
tal at low cTpw was subsequently identified in Ref. 40 

At sufficiently high pressure the most stable phase will 
be the crystal with highest density, which in this system 
is the hexagonal crystal. At sufficiently low pressure, 
the crystal structure with the lowest energy will be most 
stable, and will be the one which maximizes the patch- 
patch interactions. 

At intermediate values of the patch width, the low en- 
ergy crystals are those based on the two five-coordinate 
local environments depicted in Fig. [l] Although the 
patches are not able to point directly at those on neigh- 
bouring particles in these motifs, the loss in energy is 
not prohibitive because the angular deviations are rela- 
tively small (Fig. |3|. As the average angular deviation 
is smaller in the a local environment, the a crystal is 
slightly lower in energy than the H crystal. 

As the patches becomes narrower, the energetic 
penalty associated with the non-perfect alignment of the 
patches with the interparticle vectors increases, until a 
point is reached where it becomes favourable for the local 
environments to distort so that the patches point directly 
at three of the five neighbours. The resulting packings 
can be considered to be made of the irregular hexagon 
shown in Fig. (sjc). For the H crystal this distortion 
does not lead to any change in its space group, namely 
it remains as cmm. However, for the a crystal there is 
a symmetry breaking — there are two equivalent ways 
of dividing up the structure into these hexagons — and 
the space group changes from pAg to p2. The resulting 
crystals are virtually degenerate, with the distorted a 
structure just lower in energy. 

As the patch width is increased away from the inter- 
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FIG. 3. Deviations of the patch vectors from the interparti- 
cle vectors in the (a) a- and (b) if-type environments, (c) At 
low values of the patch width it is favourable for these struc- 
tures to distort so that three of the patches point directly at 
the neighbouring particles. The resulting structures can be 
viewed as being made up of the depicted hexagon. 



mediate values at which the a phase is stable, the ori- 
entational dependence of the potential becomes weaker, 
and there comes a point at which the hexagonal crystal 
becomes lowest in energy because of its greater number 
of neighbours. 

In the <j-phase the ratio of triangles to squares is 2. 
There are also many larger-unit cell crystal structures 
with a mixture of a and Z environments that have an 
increasing ratio of triangles to squares. We thought that 
these more complex crystals might be most stable near to 
the cr/hexagonal boundary in Fig. |2] However, we never 
found an instance where one of these structures had the 
lowest enthalpy. 



2. Annealing 
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FIG. 4. Evolution of local structural environments during 
(a) cooling and (b) annealing at T = 0.15e/c~^. p = 0.7 ea^j, 



We do not expect quasicrystalline configurations to be 
lowest in enthalpy, both because the closest quasicrys- 
talline approximants (the crystals involving both a and 
Z local environments) are never lowest in enthalpy (Sect. 
Ill A 1 ) and because the disorder associated with qua- 



sicrystals is likely to lead them to have a somewhat higher 
enthalpy than the relevant approximant. However, it 
may be that a quasicrystal is more kinetically accessible 
than possible crystals on cooling, or is even thermody- 
namically stable for a particular temperature range due 
to its greater entropy. Therefore, to search for quasicrys- 
talline behaviour we performed a series of cooling runs for 
a grid of pressure and patch width values. In these runs 
the temperature was decreased linearly from 0.5 ek~^ to 
zero over 50 000 MC cycles. Some of the resulting config- 
urations did show twelve-fold diffraction patterns char- 
acteristic of a dodecagonal quasicrystal. However, the 
peaks were very diffuse. Therefore, we subsequently an- 
nealed all the final configurations at a temperature of 
0.15 ek~^ for 10^ MC cycles. This temperature was cho- 
sen so that the mobility of the atoms was large enough to 
allow significant ordering during annealing. For the most 
part, this temperature was also below that at which the 
quasicrystals formed, the exceptions occurring near to 



the hexagonal/cr boundary in Fig. [2j where a hexagonal 
phase formed first on cooling before transforming into 
the quasicrystalline phase at lower temperature. 

Figure |4] shows the evolution of the number of particles 
with different local environments for the cooling and an- 
nealing runs for a state point that led to a dodecagonal 
quasicrystal. On cooling, the number of a environments 
increases rapidly between 0.2 and 0.1 ek~^ and is due to 
the formation of the quasicrystal. Prior to this there is 
a transient increase in the number of hexagonal parti- 
cles. Closer to the hexagonal/a boundary in Fig. [2j this 
transient increase is more pronounced because of the in- 
creased stability of the hexagonal phase. Even though 
the hexagonal crystal is not the lowest enthalpy struc- 
ture under these conditions, it is thermally stabilized by 
its large orientational entropy. 

It is noteworthy that the number of hexagonal parti- 
cles does not decrease to zero on quasicrystal formation. 
This is not due to incomplete ordering but is one of the 
features of the quasicrystalline structures. At the end of 
the cooling run, there is still a significant fraction of par- 
ticles whose local structure cannot be assigned because 
of the disorder within the configurations. 




On annealing the number of cr environments gradually 
increases at the expense of unidentified and H environ- 
ments (but not hexagonal) (Fig. [4]^b)). The resulting 
configuration is shown in Fig. [5] along with the associ- 
ated diffraction pattern. The diffraction pattern shows 
clear twelvefold symmetry and is very similar to that for a 
"perfect" dodecagonal quasicrystal produced by the "ex- 
tended Schlottmann" infiation rules^^^^^ (Supplementary 
Fig. 1), albeit with less sharp peaks. The twelvefold pat- 
tern also implies that the orientational order is coherent 
across the whole box; it could be said to be a "single 
quasicrystal" . 

The configuration is a square-triangle tiling in which 
the twelve possible orientations of the bond vectors are 
equally likely — hence, the twelvefold symmetry in the 
diffraction pattern. The radial distribution function 
shows clear peaks out to quite long range (Fig. [6]^ a)) 
because of the square-triangle order. Close inspection 
of the configuration shows that the hexagonal atoms for 
the most part do not cluster together but are instead iso- 
lated from each other and are usually at the centre of the 
dodecagonal motifs depicted in Fig. [7|^a) and (b). These 
dodecagonal motifs can join together in two ways. They 
can share an edge as in Fig. [tJ^c) (their centres are sep- 
arated by (2 + >/3)req = 4.189 (Jlj) or interpenetrate as 
in Fig. [T^e) (separation (1 + A/3)req = 3.067 (Jlj). These 
two distances are very apparent from the radial distri- 



bution function for particles in hexagonal environments 
(Fig. [6]^b)) with a clear preference for edge-sharing do- 
decagons. 

The four basic ways of locally packing the dodecagons 
that are possible using interpenetrating and edge-sharing 
dodecagons are illustrated in Fig. [7|^c)-(f), two of which 
are triangular, one of which is square and one of which 
is rectangular. These motifs can be used to construct a 
whole variety of crystal structures with large unit cells 
(and varying ratios of squares to triangles), although 
as mentioned in Sect. IIII A II we did not find an in- 
stance where these crystals had the lowest enthalpy. The 
model quasicrystal that we created using the extended 
Schlottmann infiation rules (Supplementary Fig. 1) can 
also be analysed in terms of the edge-sharing dodecagon 
motifs of Fig. [tJ^c) and (d). 

When the configuration in Fig. jsja) is examined all 
four of these motifs can be found, with the triangle 
of edge-sharing dodecagons (Fig. [tJ^c)) most common. 
When examining how these triangular, square and rect- 
angular elements begin to tile the plane, we find config- 
urations like those in Fig. [T] but now, of course, on a 
longer length scale. Furthermore, this tiling seems to be 
random and does not completely tile the plane, because 
of a significant fraction of defects, both in terms of U 
particles with an "unidentified" coordination shell and 
strings of H particles. It is noteworthy that <j and H 
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FIG. 6. Radial distribution functions for (a) all the particles 
and (b) the Z particles in the quasicrystalline configuration 
obtained at p = 0.7 ea^j, cTpw = 0.49. In (b) the peaks associ- 
ated with interpenetrating (Fig.[7]^e)) and edge-sharing (Fig. 
mc)) dodecagons are marked with arrows. 




FIG. 7. (a) and (b) Dodecagonal motifs, and (c)-(f) four 
ways of locally packing these motifs. 



crystals can form a coherent interface along the {11} di- 
rections of the a crystals, and the latter defects seem to 
often occur when the <j-like order in two adjacent regions 
is not fully in registry, and needs a line of H atoms to 
bridge the regions. 

Another source of disorder is the two possible orienta- 
tions of the central hexagon in the dodecagonal motif, as 
illustrated in Fig. [7|^a) and (b). The two forms can be 
transformed into each other by a rotation of the hexagon 
by 7r/6. For the edge-sharing dodecagons in Fig. [7| the 
orientations of the hexagons have been chosen so that 
all the particles on the edge of the dodecagons have the 
more favourable a environments. However, some relative 
orientations of the internal hexagons lead to H environ- 
ments. For example, there are two H environments at 
the shared edge if both hexagons have bonds parallel to 
the shared edge, and a few examples of such H "dimers" 
can be found in Fig. Isja). By contrast, for the interpen- 
etrating dodecagons (Fig. [rj^d) and (e)), if the hexagon 
at the centre of one of the dodecagons is rotated, the do- 
decagonal character of the other centres is lost, although 
the packing is still a square-triangle tiling. 

Another way to characterize the structures that we 
observe is in terms of the ratio of triangles to squares 
in the tiling, and the related ratio of the number of 
five- and six-coordinate atoms. For example, we find 
the ratio of triangles to squares for the configuration 
in Fig. [sja) to be 2.303, which compares with a value 
of 2 for the a crystal and 4/^3 = 2.309 for an "ex- 
tended Schlottmann" quasicrystal.'^ Similarly, the ra- 
tio of 5- to 6-coordinate environments is 13.194 for the 
configuration in Fig. [sja), which is again similar to 
the value for an "extended Schlottmann" quasicrystal, 
namely (24 + UV3)/{2 + ^3) = 12.928. These results 
clearly indicate the close similarity of the structures we 
obtain to ideal dodecagonal quasicrystals. 

So far we have only looked at the structure that results 
for one particular set of conditions. Figure [S] provides 
an overview of the structural behaviour as a function of 
pressure and patch width, and Figure [9] provides exam- 
ple final configurations for different patch widths at a 
representative pressure. 

In the top right corner of the (crpw,P) plane, i.e. larger 
patch widths and higher pressures, the hexagonal phase is 
most stable (Figure |2|. As one moves away from this cor- 
ner of the parameter space to lower pressures and patch 
widths, the onset of quasicrystal formation is signalled by 
a sharp increase in the number of a environments at the 
expense of hexagonal environments (Fig. [sja) and (b)). 
Close to this boundary, our measure of the twelvefold- 
ness of the diffraction pattern generally has high values, 
although this measure should be interpreted with caution 
as false positives can arise. For example, the superposi- 
tion of diffraction patterns from two different hexagonal 
domains leads to an anomalously high value of the twelve- 
foldness at P = 1.7 ea^j and (jpw = 0.63. The position 
of the hexagonal to quasicrystal boundary in Fig. [S] is no- 
ticeably to the left of the zero-temperature cr-hexagonal 
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FIG. 8. Dependence on pressure and patch width of the fraction of particles in (a) a, (b) Z and (c) R environments, and 
(d) \F(y — 12)1 — \F(y = 6)|, a measure of the twelvefold character of the diffraction pattern, for the final configuration after 
annealing. 



boundary in Fig. [2] because of the entropic stabilization 
of the hexagonal phase mentioned earlier. 



As one moves away from this boundary, the number 
of hexagonal atoms slowly decreases towards zero and 
the degree of twelvefoldness generally drops. Examining 
the configurations in Fig. |9] one sees a gradual crossover 
from the quasicrystalline configurations with hexagonal 
atoms at the centre of the characteristic dodecagons to 
a pure a phase in Fig.joFb). For example, about a half 
of the domains in Fig. [9[d) are pure a and the rest are 
quasicrystal-like. This crossover is driven by the hexago- 
nal environments becoming energetically increasingly un- 
favourable as the patches become narrower. 



At the narrowest patch width we considered, the low- 
energy crystals are the distorted versions of the a and 
R crystals in which three of t he patch es point directly 
at neighbouring particles (Sect. [in A 1 ). These two crys- 
tal forms are nearly degenerate, and the structure that 
results from annealing is a mixture of the two. Even in 
their distorted forms these two crystals can form coher- 
ent boundaries between them, and alternating series of 
layers of the two crystal forms can be seen in Fig. [9]^ a). 



It is also noticeable from Fig. [S] that for higher pres- 
sures at the narrowest patch width hexagonal crystals 
again form. This change reflects the energetic destabi- 
lization of the a phase at these patch widths because the 
particles can no longer form five strong interactions. 



B. 7-patch particles 

Like the 5-patch particles, the 7-patch particles have 
a local symmetry that is incompatible with crystalline 
order. However, in addition it is also physically unfea- 
sible for a particle to have seven neighbouring atoms at 
the equilibrium pair separation without particles overlap- 
ping. The question is then how do the particles manage 
to maximize their patch-patch interactions? When the 
patches are reasonably narrow, the solution the system 
finds is to use only five of the seven patches, and to adopt 
the same types of square-triangle tilings as for the 5-patch 



system. Figure 10 shows the arrangement of the particles 
in the a and R environments, and it can be again seen 
that the a environment has a smaller average deviation 
of the patch vectors from the interparticle vectors. 

Because of this tendency to form square-triangle 
tilings, the behaviour of the 7-patch system is very similar 
to the 5-patch system, and so we will much more briefly 
review the behaviour of this system. The zero temper- 
ature phase diagram shows a very similar form with a 
hexagonal crystal having lowest enthalpy at high pres- 
sure and patch width, and a a crystal at lower pressures 
and patch widths. The crossover between these forms 
occurs at slightly narrower patch widths than for the 5- 
patch system, because the 7 patches makes the potential 
somewhat closer to the isotropic limit. As for the 5-patch 
system at very narrow patch widths, this structure dis- 
torts so that three of the seven patches point directly 
at their neighbours, but this time the acute angle in the 
hexagonal units is 51.43° not 72°. 
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FIG. 10. Deviations of the patch vectors from the interpar- 
ticle vectors in the (a) a- and (b) iif-type environments for 
7-patch particles. 



Again close to the hexagonal/<j boundary our anneal- 
ing simulations result in quasicrystalline configurations 
with clear 12-fold diffraction patterns (Fig. 11). Inspec- 
tion of the configurations also again shows the same do- 
decagonal motifs (Fig.[7|) as found for the 5-patch system. 



IV. CONCLUSIONS 

Here we have presented results for a simple two- 
dimensional patchy particle system that exhibits a rich 
ordering behaviour, in particular showing the formation 
of dodecagonal quasicrystals based on square-triangle 
tilings for certain parameter ranges for particles with 
five and seven regularly arranged patches. These ex- 
amples can be added to the increasing list of model 
systems that have been shown to form quasicrystals in 
simulations.^^lf^ The system also provides a nice model 
system that illustrates how the local structural propen- 
sities can lead to complex globally-ordered structures. 
In particular, quasicrystals were observed in the region 
of parameter space near to where the enthalpically pre- 
ferred structure changes from five-fold to six-fold coordi- 
nation, because of the presence of both coordination en- 
vironments in the quasicrystal. Interestingly, the struc- 
ture of the quasicrystals is quite similar to model square- 
triangle dodecagonal quasicrystals produced by inflation 
rulesPT^^ 

One interesting question concerning the quasicrystals 
that we observe is whether they are thermodynamically 
stable or just a kinetic product. As the enthalpy dif- 
ference between the quasicrystal and the crystal is not 
that large, it is not unreasonable that the quasicrystal 
could become thermodynamically stable due to the en- 
tropy associated with the many possible configurations 
for the quasicrystal. However, to confirm this hypothesis 
would require the free energy of the quasicrystal (or the 
free energy difference between it and the crystal) to be 
computed, but it is not clear how this could be achieved. 
During some of the cooling simulations, the system passes 
from a liquid to a hexagonal crystal and then to a qua- 
sicrystal as the temperature is decreased, showing that 
there is a parameter range where the quasicrystal is more 



stable than the hexagonal crystal. However, on heating 
the a crystal, we have never seen it transform into a qua- 
sicrystal, but instead it directly melts. 

An interesting contrast to the behaviour we see here 
is provided by particles th at ha ve five-fold symmetry in 
their repulsive int er act ions ,1^^^^ rather than their attrac- 
tions. Although hard pentagons show an interesting set 
of solid phases structures, none of them are quasicrys- 
talline. 

Another important question is whether a system that 
shows behaviour similar to our model could be experi- 
mentally realized. Methods to synthesize patchy colloids 
and nanoparticles are developing rapidly, but as far as 
we are aware, there are none yet available that could 
produce, say, particles with a five- fold symmetric dis- 
tribution of patches. By contrast, effectively restricting 
such colloidal systems to two-dimensions would be rela- 
tively straightforward. For example, if there is a density 
mismatch between the colloids and the solvent, sedimen- 
tation can lead to monolayer formation at the base of the 
sample, as has been done in recent experiments where a 
two-dimensional Kagome lattice was assembled from 'tri- 
block Janus' particles.'^ 

Another possible system in which these dodecagonal 
quasicrystals could be potentially realized is the multi- 
arm DNA motifs produced by the group of Chengde 
Mao .I^^H^ Each arm is made of two-parallel double helices 
and both have dangling single-stranded ends that allow 
them to bind to other such motifs with a well-defined rel- 
ative orientation. These effective torsional constraints on 
the interactions lead to quasi-two-dimensional growth be 
it into sheet s^^ ^^ or closed shells. Interestingly, simi- 
lar to the equivalent patchy particles, the 3-^ 4- and 6 -arm 
motifs form two-dimensional honeycomb j ^^ * ^^ * squar j^^^^ 
and hexagonaf^ lattices. Furthermore, the five-arm mo- 
tifs can form a a-phase like lattice.^^ 

However, there are also a number of significant dif- 
ferences between the DNA multi-arm motifs and patchy 
particle systems. Firstly, the central loop to which the 
stiff double helices are connected has a certain degree of 
flexibility allowing the relative angles of the arms to vary. 
By contrast, the positions of the patches in our model are 
fixed. 

Secondly, the "valence" of the DNA motifs is fixed; 
i.e. the five-arm motifs can only ever bond to five other 
such motifs. By contrast, for our patchy particles at in- 
termediate values of the patch width, the particles can 
adopt five- fold or six- fold coordination environments; it 
is this flexibility that allows the system to form the qua- 
sicrystals. Therefore, if the DNA motifs are to be able 
to form a quasicrystal, a mixture of 5- and 6-arm motifs 
with the right stoichiometry would be required. For an 
ideal dodecagonal quasicrystal, the required ratio of 5- 
and 6-arm tiles would be 12.928 : 1. However, even with 
such a mixture, it might be that the system prefers to 
phase separate into a and hexagonal crystals, similar to 
what has been seen for mixtures of 3- and 4-arm tiles.l^ 

In future work, we tend to explore this possibility 
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FIG. 11. (a) Configuration after annealing and (b) and (c) associated diffraction patterns (for two different ranges of q) at 



p — 0.7 ecr^ 



0.42 for a system of 7-patch particles. 



further, first by using a coars e-gra ined model of DNA 
that we have recently developecP^l^to characterize these 
DNA multi-arm motifs, including their flexibility and the 
angular specificity of their interactions. Secondly, this 
information will then be used to create a patchy-particle 
representation, where the patch positions are not rigidly 
fixed but are constrained by an internal potential. 
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V. SUPPLEMENTARY INFORMATION 

A. Model quasicrystal 

As a comparison to the diffraction patterns we have 
calculated from our simulated configurations, here we 
report the diffraction pattern from a model dodecago- 
nal quasicrystal. Dodecagonal tilings can be created by 
projection from a four-dimensional space but usually in- 
clude rhombs as well as squares and triangles. An al- 
ternative approach to create dodecagonal square-triangle 
tilings is to use an infiation approach where an increas- 
ingly large structure is iteratively built up by apply- 
ing rules for scaling simpler motifs (in our case squares 
and triangles) and then replacing them by more com- 
plex ones.^^ 66 68 HQYe we use what has been termed the 
"extended Schlottmann" infiation rules .1^^^^ The effect of 
these rules is at each iteration to replace each vertex in 
the square-triangle tiling by one of the dodecagons in Fig. 
[7|^a) and (b) depending upon the local geometry of the 
vertex. Supplementary Figure 1(a) shows the structure 
generated after 3 iterations of these rules starting from 
an initial centred hexagon. The corresponding diffraction 
pattern has clear 12-fold symmetry (Supplementary Fig. 
1(b)) and a very similar pattern of peaks as those in Fig. 
|5]and[lll 
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